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We study mechanism of magnetic energy transport, motivated by recent measurements of the 
thermal conductivity in low dimensional quantum magnets. We point out a possible mechanism 
of enhancement of the thermal conductivity in gapped magnetic system, where the magnetic 
energy transport plays a crucial role. This mechanism gives an interpretation for the recent 
experiment of CuGeOa, where the thermal conductivity depends on the crystal direction. 
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§1. Introduction 

The relationship of the thermal conductivity and magnetic state in various low dimensional quan- 
tum magnets has been attracted interests.i'i'il&i'& The thermal conductivity can be a strong 
experimental instrument to investigate the magnetic properties in materials. For the magnetic sys- 
tems with large exchange coupling constants, the magnetic energy transport dominantly contributes 
to the thermal conductivity. Low-dimensional spin systems with many conserved quantities tend 
to possess a ballistic energy transport, which means an infinite thermal conductivity as studied in 
the theoretical side.i'i'& These facts can give a possible mechanism of a large magnetic thermal 
conductivity. Solgubenko et al. observed such nondiffusive magnetic energy transport in SrCu02 
and Sr2Cu03 which are described by the quasi-one dimensional isotropic antiferromagnetic Heisen- 
berg model.0) These materials have very large exchange coupling constant J ~ 2000i^' so that the 
magnon runs through the system even at 100 K. They found that the energy is transmitted by 
magnon at this temperature. 

Unusual enhancement of conductivity is observed at low temperatures below a spin gap temper- 
ature in gapped spin system such as the two-dimensional dimer spin system SrCu2(B03)2, and the 
ladder system (Sr, Ca)i4Cu2404i.i'i'&i'il' This enhancement is attributed to the contribution of 
phononic energy transport. This phononic energy transport is quantitatively explained in view of 
'resonant scattering process' between phonon and spin gap.Bi'E^ That is, the magnetic system is 



excited absorbing a phononic energy, and another phonon of the same energy is emitted deexciting 
the magnetic system at the same time.0^ These experiments suggest that at high temperatures the 
magnetic energy transport dominantly convey the heat, while at low temperatures the magnetic 
gap suppresses the magnetic energy transport. Thus phononic energy transport become dominant. 
Therefore the energy carrier changes from spins to phonons around the spin gap temperature. 

The spin-Peierls (SP) system is another gapped system where the energy gap is formed below a 
critical temperature Tgp. CuGeOs is a typical material which shows the SP transition. Ando et 
al. studied the thermal conductivity in the chain direction in CuGeOs and found the enhancement 
of the conductivity below Tgp0) which is similar to that found in the two-dimensional dimer and 
ladder systems.i'i'0i'& Ando et al. attributed this unusual enhancement of conductivity to the 
increase of the phononic transport due to the decrease of spin-phonon coupling interaction. The 
similar enhancement is also observed more clearly in NaV2 05 where the charge ordering transition 
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occurs J 

Recently Hofmann et al. compared the properties of thermal conductivity of SrCu2 (603)2 and 
CuGeOs.© They investigated the thermal conductivity along different crystal axes, although the 
previous experiment in CuGeOs was done only in the chain direction.0) Generally temperature 
dependence of phononic thermal conductivity does not depend on the crystal direction, whereas a 
magnetic one does on it. Thus by studying the crystal direction dependence of thermal conductivity, 
whether the enhancement of the thermal conductivity is caused by phonons or not is clarified. 
SrCu2(B03)2 is described by the Shastry-Sutherland model0) which has the exact dimerized ground 
state and almost dispersionless triplet excitations where the magnetization is localized. Thus below 
the spin gap temperature, we expect that the phononic energy transport is dominant because of 
small magnetic contribution. Actually this prospect is supported by the observation that the 
thermal conductivity does not depend on the crystal direction. On the other hand, they found that 
CuGeOs shows the direction dependence of thermal conductivity. 

In this paper, motivated by these recent experiments, we theoretically investigate the role of 
magnetic transport when the system has the energy gap, and propose a mechanism of enhancement 
of the thermal conductivity in magnetic systems. We consider a bond-alternating antiferomagnetc 
Heisenberg chain, and investigate the thermal conductivity focusing on the dependence on the 
energy gap. By the quantum master equation approach, we study a spin system which is connected 
to two reservoirs of different temperatures. We found that the thermal conductivity is enhanced 
due to the presence of bond-alternation around the spin gap temperature. This observation gives 
an interpretation that a small bond-alternation yields an enhancement of specific heat but does 
not drastically change the mean free path. This interpretation is consistent with Hofmann's recent 
experiments. 

This paper is organized as follows. In §2, we introduce the model, and the method to investigate 
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the thermal conduction is explained in §3. Summary and discussions are given in §4. 
§2. Model 

The system we shall consider is an alternate Heisenberg spin chain described by, 

= Jj] (l-(-lf,5)s,-S,+i, (2.1) 

£=l 

where J is an exchange interaction which is taken as a unity throughout this paper, and is 
the ith spin. The number of spins N is taken to be even. A finite value of parameter 5 causes a 
bond-alternation, which yields the dimerized ground state and a finite spin gap0 We study the 
thermal conductivity for various values of 5. 

The thermal conductivity n{T) at the temperature T is generally given by the Green-Kubo 
formula which reads as, 

where {., .} means the anti-commutation relation, and (...) is the equilibrium average at the tem- 
perature T. The operator J{u) is the total heat current operator at the time u in the Heisenberg 
picture, i.e., J{u) = exp{iTCu)J ex.p{—i7iu). The total current operator J for (2.1) is calculated by 
the continuity equation of the energy as, 

Af-2 

J = -iJ ^ (l- (5^) [S^ • Sf+i, S£+i • S^+2] 

N-2 

= -j(l-6^) J2 {Sex Se+i) -8^+2 (2.3) 
e=i 

where [., .] means the commutation relation. In the case of the isotropic Heisenberg model with the 
periodic boundary condition, i.e. 5 = 0, the total current (2.3) is a conserved quantity. Therefore 
the Green-Kubo formula (2.2) diverges, which means that magnons are not scattered and the mean 
free path of each magnon diverges. In real magnetic materials, the relaxation time becomes finite 
due to the other sources such as spin-phonon scattering and impurities, and so on. 

In order to take these sources into account for the isotropic case 5 = 0, more practical formulation 
of the thermal conductivity is used especially when the magnetic low excitation contributes to the 
energy transport at low temperaturesjillli^ 



k{T) = vI I dkV{k)T{k)^ 



huj{k) 



(2.4) 



where uj{k) is the low energy dispersion which is taken to be the Des Cloizeaux-Pearson mode0^ 
uj{k) = ^ sin k. vq is the group velocity. 'D{k) and T{k) are the density and the relaxation 
time of the mode k, respectively. Eq. (2.4) can be also written as k = J2k^kiT)Dk{T) where 
Ck{T) = T>{k)-^ ^i3H'^(k)_^_i aiid Dk{T) = VQT{k). Ck{T) and Dk{T) are interpreted as the specific 
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heat and the diffusion constant of energy for the mode k, respectively. This is derived with the same 
spirit in the phononic transport discussed by the Peierls-Boltzmann equation.© The relaxation 
time T{k) is determined by various origins such as impurities and spin-phonon coupling, and so on. 
Practically, the relaxation rate is calculated by the Carraway's method,© i.e., by the summation 
over relaxation rates originated from these origins. In this formula, the mean free path Ik is defined 
using the relaxation time t(A:) as Ik = VQT{k). In CuGeOs, above Tsp where the Hamiltonian is 
roughly regarded as the isotropic Heisenberg chain 5 = 0, the mean free path is estimated as 500 
lattice size. We should note that the mean free path is a dynamical property and is different 
from the spin correlation length in the equilibrium state. 

In order to apply the formula (2.4) to the case of finite (5, we must know the low energy dispersion 
with finite energy gap, although we do not know exact analytical expression for it. A finite energy 
gap causes the localization of the magnetic excitation which suppresses the energy diffusion, whereas 
energy gap enhances the heat capacity around spin gap temperature. For instance, when (5 = 1, 
the system is completely separated into N/2 local spin pairs, and the group velocity is 0. Therefore 
this extreme case gives zero thermal conductivity due to vanishing mean free path. Thus it is not 
trivial whether the existence of energy gap enhances thermal conductivity or not. In this context, 
we will study the thermal conduction by a direct numerical method. 

§3. Method and Results 

We investigate how the thermal conductivity behaves as a function of 5. It is difficult to treat 
the Green-Kubo formula (2.2) numerically because it requires infinite N . Therefore the quantum 
master equation approach is more tractable in numerical simulations. We study the spin system 
which is directly connected to the reservoirs of different temperatures at the ends. We here consider 
a system of = 6 and N = 8. Although this lattice sizes are small, we believe that the essential 
mechanism of magnetic energy transport is clarified. The time-evolution of the system is determined 
by the following master equation (e.g., see the reference and references therein); 

^ = -i [n, p{t)] - A {CLPit) + CKpit)) , (3.1) 

where the first term in the right-hand side corresponds to the pure quantum dynamics of the 
system, and >Cl and jCr, express the contact with the thermal reservoirs of the inverse temperature 
/3l(= I/Tl) at the left end and /3r(= 1/7r) at the right end, respectively. The parameter A is the 
coupling strength. The dissipative term (a = L, R) is given by 

CMt) = {[Xa, RaP{t)] + [A„, ii„p(t)]t) , (3.2) 

where Al and Ar are the system's operators directly attached to the left and right reser- 
voir, respectively. Here we take Al = Sf and Ar = Sf^. The operator Ra is given by 
{k\Ra\m) = {Ek — Em){e^°'^^''~^'^^ — l)^^(fc|AQ,|m) in the representation diagonalizing the Hamil- 
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tonian (2.1) as T-L\k) = Ei^\k) and T-L\m) = Em\m). If /?l = /3r, the master equation (3J) gives 
the canonical distribution as the stationary solution, which is easily checked by substituting the 
canonical distribution for the density matrix p{t). We numerically integrate the equation (|3.1D and 
obtain the stationary density matrix pst starting from the initial density matrix of the canonical 
distribution at temperature Tr. The simulation was carried out by the fourth order Runge-Kutta 
method with the time step 0.01 with A = 0.01. 
We define the thermal conductivity by measuring the stationary energy current Tr [pst'Pj setting 

Tr = Tl + AT, (3.3) 

where AT is a temperature difference between two reservoirs, and is taken as AT = 0.3. Thereby 
the thermal conductivity k(T) is defined as 

Tr(jp, 



k(T = Tl) = ^^. (3.4) 

We present the results of the thermal current k(T) AT normalized by N — 2, which is the number 
of terms in the summation (2.3) In Fig. 1(a) and (b), the normalized thermal currents are shown 
for = 6 and = 8, respectively. In the inset, the data of the specific heat are also shown. 
The normalized thermal current show little quantitative difference between the cases of = 6 and 
N = 8. In the figures, we obtained finite thermal conductivities even in the isotropic case, 6 = 0. 
As explained in previous section, it must diverge in the Green-Kubo formula. These finite values of 
thermal conductivity are attributed to the finite size effect in the presence of the thermal reservoirs 
at ends of the system, and the conductivities must diverge in the limit of A^ — > oo. Actually the 
conductivity for A^ = 8 is larger than that for A^ = 6. This size dependence is also observed for the 
cases of finite values of 6. 

In the isotropic case, 6 = 0, the overall form of the thermal conductivity is similar to that obtained 
by Kliimper and Sakai who exactly calculated the amplitude of zero frequency contribution in the 
Green-Kubo formula@ That is, the thermal conductivity has one peak at about T ~ 0.5. This 
behavior is a common characteristic which does not change when the system size increases. At very 
low temperatures k(T) should be in proportional to T due to the temperature dependence of specific 
heat, which is similar to the Casimir's theory in phononic transport.!!) This low temperature 
property, however, cannot be checked for such finite sizes. We must note that even in the isotropic 
case 6 = 0, the system has a finite energy gap due to the finite size effect. When 6 becomes finite, the 
gap increases, e.g., in the case of A^ = 8, A^; = 0.39269, A^; = 0.74750, 1.09041, 1.41106, 1.71290, 
and 1.85807 for 6 = 0.0,0.2,0.4,0.6,0.8, and 0.9, respectively. 

Small bond alternations enhance the thermal conductivity in spite that the alternation tends to 
separate the system into local spin pairs. This enhancement is observed for even 6 = 0.8, and this 
feature is quite robust. In the case of very large 6, the thermal conductivity is reduced due to the 
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effect of separation of the system. The temperature of the peak of thermal conduction roughly 
corresponds to half of energy gap like the Shotky-type specific heat. Actually the overall behavior 
of thermal conductivity is similar to that of the specific heat. Thus from this observation, we 
conclude that the mean free path does not drastically changed by the bond-alternation. We expect 
that the unusual enhancement of the thermal conductivity in the spin-Peierls system is caused by 
magnetic energy transport when Tsp is the same order of A£?. That is, in the spin-Peierls system, 
k{T) follows the data of (5 = at the high temparture region, whereas below Tsp, S suddenly 
increases and k{T) is given by k,{T,5 ^ 0), which gives a sharp peak when Tsp is the same order 
of AE. In CuGeOs, Tsp is the same order of AE, i.e., Tsp ~ UK and AE ~ 2QK. 

§4. Discussion 

Magnetic energy transport is becoming a strong instrument to capture the magnetic properties 
of quantum magnets. In order to realize an accurate interpretation for obtained experimental data, 
theoretical arguments are necessary. We demonstrate how the thermal conductivity is enhanced 
due to the magnetic interaction, and proposed that besides the enhancement due to the increase 
of phononic energy transport, magnetic energy transport also causes an enhancement at a tem- 
perature of the order of /S.E. This enhancement is directly related to that of the specific heat. 
This correspondence indicates that the mean free path does not drastically changed by a small 
bond-alternation. If the bond-alternation disturbs the dynamics of the system and the mean free 
path becomes much smaller than the isotropic chain, this correspondence between the thermal 
conductivity and specific heat is no more valid. The numerical simulation, however, indicates that 
there is some region where this mechanism of enhancement really exists. We hope that this effect 
by magnetic energy transport is also studied in experiment. 

So far, it is difficult to control the heat flow in material. However making use of the magnetic 
energy transport, we may plan a new device that controls the energy flow. In such device, the 
energy flow could be drastically changed by controlling some parameters such as external magnetic 
fields. The creation of a energy gap provides a promising mechanism which enhances the energy 
flow. It is interesting to study mechanisms to control magnetic transport in more detail, which will 
be reported elsewhere. 
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Fig. 1. The thermal conductivities for various 6 with the system size, (a): A'^ = 6 and (b): A'' = 8. The inset is the 
specific heat. 
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